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I. INTRODUCTION 


The main motivation of this work is to understand, in rigorous terms, the Hamiltonian 
formulation and subsequent Fock quantization of linear systems consisting of fields coupled 
to point masses (and, eventually, other low-dimensional objects). An important question 
that we will answer in the paper regards the possibility of writing the Hilbert space of such 
a compound system as a tensor product of Hilbert spaces associated with the point masses 
and the field. The point masses are introduced to model external devices that can be used 
both to excite the system and to act as detectors sensitive to the “field quanta”. In this last 
sense they can be thought of as generalizations of the Unruh-DeWitt particle detectors and 
similar devices, used in the discussion of quantum field theories in curved space-times and 
accelerated frames [3-0|. We wish to emphasize from the start that the point masses -that 
we introduce already at the classical level- have nothing to do, in principle, with the field 
quanta despite the fact that the latter are usually interpreted as particles or quasiparticles. 


The types of systems that we consider here are related -but not equal- to field theories 
defined in bounded spatial regions and share some features with them, as will be explained 
in the following. There are several important field theoretic models that display interesting 
and non-trivial behaviors when defined in such regions (or more generally in the presence of 
spatial boundaries). Among them we would like to mention Chern-Simons and Maxwell- 
Chern-Simons models in 2+1 dimensions |HM3] , Yang-Mills theories in 3+1 dimensions |8;], 
and general relativistic models such as the isolated horizons used to study black hole entropy 
in loop quantum gravity jo|. 10]. Part of our work is motivated by the comments appearing 
in EH regarding the use of boundaries with the same purpose as the classical point-particles 
that we introduce here. 


The standard approach to derive the Hamiltonian formulation for field theories, especially 
when gauge symmetries are present, relies on the methods developed by Dirac [12|. These 
are straightforward to use in the case of mechanical systems with a finite number of degrees 
of freedom and -in simple circumstances- can be adapted to field theories if one is willing 
to accept a certain lack of mathematical precision (regarding, for instance, the functional 
spaces describing the field degrees of freedom). In the presence of boundaries, however, the 
naive implementation of the Dirac algorithm is awkward and often leads to incomplete or 
plainly wrong results (see jl3j and references therein). This is even more so for the model 
that we consider in the paper where we have both boundaries and particle-like objects 
associated with them. Some technical details regarding the difficulties of implementing the 
Dirac algorithm for field theories can be found in [l3 |. 


These difficulties -and other important ones- can be avoided by the use of the geometric 
constraint algorithm developed by Gotay, Nester and Hinds TH~16j- This method provides 
a rigorous, geometric and global way to obtain the Hamiltonian description of field theories 
and pays due attention to functional analytic issues. In particular it provides a completely 
detailed description of the spaces where the Hamiltonian dynamics takes place and can thus 
be used as the starting point for quantization. This is especially useful for linear theories for 
which Fock quantization can be rigorously defined starting from the complexification of an 
appropriately defined real Hilbert space of classical solutions (along the lines described, for 
example, in jl7j hi the context of quantum field theory in curved spacetimes). The precise 
construction of the Fock space is important in order to discuss the eventual factorization of 
the Hilbert space V. of the system in the form V. = "Hmasses <8> 'Hfieid that would account for 
a clean separation between quantum point particle and field degrees of freedom. 
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In the present paper we will study a model consisting of a finite length elastic string 
attached, at the ends, to point particles subject to harmonic restoring forces (in addition 
to the ones exerted by the string). Although similar systems have been considered in the 


literature [18|, [19[ our approach will concentrate on several mathematical issues —relevant 
from the physical point of view— that have not been discussed elsewhere. Specifically we 
will introduce a natural way to deal with this model by using a particular class of measures 
and the Radon-Nikodym (RN) derivatives defined with their help. This measure theoretic 
approach is suggested by Mardsen and Hughes in [ 20 ], page 85]. 

The methods and ideas put forward here can be exported to more complicated situations. 
For instance, it should be straightforward to generalize them to deal with higher dimensional 
systems where, in addition to point particles, other low-dimensional objects could be coupled. 
From the perspective of the standard quantum held theory in curved spacetimes, we are 
extending the usual approach by replacing the all-important Laplace-Beltrami operator by 
an elliptic operator defined with the help of a measure that combines physical features of 
both the held and point masses of the system. We would like to remark here that, as free 
theories are essential building blocks for the perturbative quantization of non-linear models 
(for instance, they play a central role to define the Fock spaces used in their description), 
it is important to understand them well as a hrst step to consider their quantization in the 
presence of boundaries and/or lower-dimensional objects. 

The layout of the paper is the following: After this introduction, section |H] will be devoted 
to the classical description of the system under consideration. In particular, in lll.2l we present 
a short discussion of the resolution of the evolution equations by separation of variables. 
This motivates the introduction of non trivial measures that play a relevant role in section 
111.31 where we provide an alternative Lagrangian formulation. In section III. 41 we obtain 
the Hamiltonian description by implementing the GNH algorithm [14-16]. We provide the 
precise description of the functional spaces relevant for the model. As it will be shown 
these are generalizations of Sobolev spaces that can be understood in a neat way with the 
help of a scalar product defined in terms of appropriate measures. I 11 section [TTT] we use 
the Hamiltonian description of the system to build a Fock space and quantize. We will pay 
particular attention to the characterization of this Fock space as a tensor product of Hilbert 
spaces associated with the held and the particles. We end the paper with our conclusions 
in section [TV] and one appendix where we give a number of useful mathematical results. 

The reader interested only in the physical results of the paper should skip sections III. 31 
and 111.41 where a number of mathematical details are provided, and go directly from section 
III. 21 to the quantum description in terms of physical modes discussed in section Hill 


II. CLASSICAL DESCRIPTION 
II. 1. Lagrangian and field equations 

Let us consider a model consisting of an elastic string of finite length (in 1 + 1 dimensions) 
coupled to two point masses located at the ends and attached to springs of zero rest length. 
Both the string and the masses are subject to restoring forces proportional to the deviations 
from their equilibrium configurations. For definiteness we will consider that the motion of 
the system is longitudinal, although this is not essential. From a logical perspective the 
equations of motion for such a system should be obtained by analyzing with due care the 
forces acting on the string and the masses, and using Newton’s laws. In practice, however, 
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it is more convenient to use an action written in terms of an easily interpretable Lagrangian 
and formally derive them by computing its first variations. Let us, then, start from the 
action ^ 

S(u) = I L(u(t),ii(t)) dt (II.1) 

Jti 

where the Lagrangian, for smooth Q and V, has the following form 


LIQ.V) 



( 11 * 2 ) 


In the preceding expression (•,•) denotes the usual scalar product in the Hilbert space L 2 (0, £) 
defined with the help of the Lebesgue measure /i^, £ is the length of the unstretched string, 
A its longitudinal mass density, e its Young modulus, m 2 > 0 is the spring constant per unit 
length associated with the restoring force acting directly on the string, M$ and Mi are the 
masses of the point particles and ko , ki the elastic constants of the springs attached to them. 
The field Q(x) represents the deviation of the string point labelled by x from its equilibrium 
position. Spatial derivatives are denoted by primes and time derivatives by dots. 

As it can be seen the Lagrangian has terms of “field” and “particle” types involving 
spatial derivatives of first order, at most. It is convenient at this point to choose units 
of length, time and mass such that £ = e = A = 1 (which in particular implies that the 
speed of sound is c 2 := e/A = 1). Notice that by doing this we exhaust all the freedom 
in the choice of units so we will have to keep h explicit when quantizing the model. To 
remind the reader of this choice we will rename the remaining constants in the Lagrangian 
as uj 2 := m 2 /e,tj 2 := kj/Mj , /ij := Mj/e so that (III. 2D becomes 


L(Q,V) = \(V,V) - - y(Q,Q) 


E f(w)’-W 


je{0,i} 


(IL3) 


It is interesting to mention at this point that the positions of the point particles are not 
classical independent degrees of freedom as they are given by continuity by the position of 
the ends of the string. This fact will have an analogue in the quantum description of the 
model (non-factorization of the Fock space, see section uni. 

The equations of motion derived from (III. 3D are: 

u(x, t) — u"(x, t) + uj 2 u(x, t) = 0 , 16 (0,1), (II.4) 

/x 0 w( 0 , t) — u'{ 0, t ) + [i o uJqu(0, t) — 0 , (II.5) 

/ ui'u(l, t) + u'{ 1, t) + n , t) — 0 . (II. 6 ) 


As can be readily seen the time evolution of the deformation of the string is governed by the 
1+1 Klcin-Gordon equation whereas the point masses at the boundary points move under 
the combined force exerted by the springs and the string (given by the spatial derivatives at 
the boundary). It is important to notice at this point that the preceding equations are not 
conventional in the sense that (1II.5[) and (III .61) are not standard boundary conditions because 
they involve second order time derivatives. We will see in section III. 21 that this feature 


1 We follow the custom of calling this constant m 2 because the squared mass of the quantum excitations of 
the usual Klein-Gordon field is proportional to it. The case m = 0 differs slightly from the one that we 
study but can be approached with the same methods. 
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qualitatively changes the type of eigenvalue problem that has to be solved to identify the 
normal modes and characteristic frequencies. Actually it renders the present problem quite 
non-trivial because the relevant eigenvalue equations are not of the standard Sturm-Liouville 
type and hence the classical theorems that ensure that the normal modes form a complete 
set cannot be applied. There is however a workaround to prove such result that consists in 
introducing a different measure space. We anticipate now the result of this approach: 

Proposition 1. The equations (III.4D - flII.6D are contained in the equations: 

u — (A m — uj 2 )u = 0 x G [0,1 ], p-a.e. (II.7) 

■ c\n 

(-l)^(j) - A(j)u(j) = 0 je {0,1} , (II.8) 

that consist of a 1+1 dimensional Klein-Gordon equation on the interval [0,1] subject to 
Robin boundary conditions written in terms of the Radon-Nikodym derivative A- with respect 
to the measure p = a^Q+pL+aidi. The parameters a 3 are related to the physical parameters 
of the problem by aj{l — ajPjiuj 2 — uj 2 )) 2 = p 3 , A(j ) := (uj 2 — uj 2 )^/pja 3 , and the Laplace-like 
operator 

A " :=(1+C) c^5 

is defined in terms of the function C(j) := A(j)a 3 , C(x ) := 0 if x ^ 0,1. 

In the previous proposition p-a.e. stands for “/i-almost everywhere”, (i.e. the equality 
can fail to be true in a set of zero p-measure, at most). The measure p is defined in the 
Appendix (see Equation IA.2|) and the domain of A^, which specifies the regularity conditions 
on the solutions to the equations of motion, is given in Appendix IA.III1 Instead of giving 
a direct proof of this result, we will apply separation of variables to (III.4D - (III.6D . In this 
process some issues will arise, the most important one being that the standard Laplacian is 
not self adjoint. This will lead to the introduction of a new self-adjoint Laplace-like operator 
(in an appropriate functional space) in terms of which we obtain equations (III.7D and (III.8D . 


II.2. Classical description of the model: solving the field equations 

In this section we consider the resolution of the equations of motion (1II.4D - (1II.6D by using 
the method of separation of variables. By writing u(x,t) = X(x)T(t) we get 


f = (A-<h 2 )T 

(11,9) 

X" = XX 

(II. 10) 

X'(0) = p o (X+u 2 o -u 2 )X(0) 

(11.11) 

X'(l) — —pi(X + u\ — o; 2 )A(1) 

(11.12) 


where AgR. In this form these equations do define an eigenvalue problem for A" with one 
key (and relatively unusual) feature: the eigenvalue appears also in the boundary equations 
(III. IIP . (III. 1211 . This means that we are not directly dealing with a Sturm-Liouville problem 
and, hence, we cannot directly import the usual results that characterize the eigenvalues A 
(do they exist? are they isolated? are they bounded?) and the corresponding eigenfunctions 
(are they a complete set? are they orthogonal?). The answer to these questions is important 
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in order to expand the general solution to the equations of motion as a functional series of 
eigenfunctions and also to quantize the system. 

In any case, a lot of information can be gathered in practice by solving the concrete 
eigenvalue problem that we have at hand so we sketch now the computation of the eigenvalues 
and the eigenfunctions. 


i) Negative eigenvalues A = —uj 2 < 0 

The solutions to (III. 101) are of the form X\ = acos(cnx) + 6sin(a;x) with a, b G R. The 
conditions (III.111) . (III. 121) imply that 

/a 0 (uj 2 — Acufya + uib = 0 

[jiiiuj 2 — Au3 2 ) cos a; + uj sin uj ^ a + 

3 we have int: 
for a, b if and only if 


UJ 


Acn 2 ) sin oj — uj cos uj j b = 0 


Where we have introduced the shorthand Acn 2 := u 2 — u 2 . These have non-trivial solutions 


[uj 2 — p 0 Pi(w 2 — AuJq)(uj 2 — Acn 2 ) j sincn+ [hq(uj 2 — AcUq) + fii(uj 2 — Aca 2 ) j uj coscn = 0. (11.13) 


It is straightforward to see that (III. 131 ) has an infinite number of solutions for uj for every 
(physical) choice of parameters /ao , Hi ,uJq ,ujf ,UJ 2 ] in fact, there exists no G N, such that 
every interval of the form (kir, (k + l)7r) contains one and only one solution of (III. 131 ) for 
every fceN, k > no ■ In the asymptotic limit k —» oo we have 


uj k = kn + 



(11.14) 


so we see that we actually have an infinite set of negative eigenvalues and also that they 
grow without bound. Finally the eigenfunctions have the form (labelling them with k and 
with a minus superscript to indicate that the eigenvalue is negative) 

X^(x) = uj k cos (uj k x) + p 0 (Acho - ujI) sin(cufccr). 


ii) Positive eigenvalues A = uj 2 > 0 

In this case the eigenfunctions X\{x) have the form X\(x) = ae + be~ ux where the real 
coefficients a, b must satisfy now the conditions 

— /io(cn^ + Acug) j a — ^cn + po(^G + Acuq)^ 6 = 0 

[uj + pi(cn 2 + Ach 2 ) j a — e~ u [uj — pi(cu 2 + Ach 2 ) j 6 = 0 

that have non-trivial solutions if and only if 

e w [uj — /^o [uj 2 + Acn5)j [uj — /ii(ca 2 + Acnf)^ 

—[uj + Hq(uj 2 + Acug) j [uj + pi(cn 2 + Acn 2 )^ = 0 . 

It can be seen that this equation has a finite number of solutions N (maybe none) depending 
on the particular choices of the physical parameters defining the problem (e.g. if both 
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A uij > 0, there is no positive eigenvalue). Notice that, as the energy is constant, the 
function T cannot have an exponential growth so a < uj 2 and, in fact, the limit case a = uj 2 
(where T has a linear behavior) happens only if ujj — 0 = uj (in the paper we are assuming 
uj > 0). Consider now k G {—1,..., — N} as a negative and finite counter, the corresponding 
eigenfunctions are 

X£(x) = (u k + /j. 0 (ul + Aul))e UkX + (c o k - p 0 (u 2 k + Auj^))e~ UkX . 


iii) Zero mode A = 0 

It is easy to check that A = 0 is an eigenvalue if and only if (1 + /i 0 Au; 2 )(l + piAcof) = 1, 
in which case the solution is simply X°(x) = A°(0)(1 + HqAuJqx). This can only happen if 
AuJq = 0 = A uj\ or AcDq < 0 < Auj\ or Auf < 0 < Au3q. 

From now on we will collectively denote the eigenvectors as X n with X n = A+ when 
n G { — 1,..., — N}, X 0 = X° (“when it exists”), and X n = X~ when n G N. Their 
associated eigenvalues A„ are, respectively A n = uj ^ when n G { — 1,..., — N}, A 0 = 0, and 
\ n = —uj\ when n G N. 

Notice that, at this point, we have found the solutions of the form u n (x,t ) = X n (x)T n (t) 
to equations (1II.9D - (III. 121) . In order to find all the solutions to (III. dD - (III. 611 we have to prove 
that the eigenfunctions {A^} form a complete set in an appropriate functional space. This is 
not straightforward because, as mentioned before, we do not have a standard Sturm-Liouville 
problem. 


The generalized Sturm-Liouville problems of the form defined by (III. 10D - (III. 12f) have a 
long history both in mechanics and mathematics (see, for instance [181, 21h23I] and references 
therein). A direct but important observation is the fact that eigenfunctions X k , correspond¬ 
ing to different eigenvalues A*,, are not orthogonal with respect to the standard scalar product 
in L 2 (0,1) but are orthogonal with respect to the following modified scalar product 


{(u, v)) /j o u(0)v( 0) + piu{l)v{l) + / u-vp L . 


(11.15) 


This can be readily proved by taking two such eigenfunctions X m , X n (associated with 
different eigenvalues A m , A„), integrating the following identity over the interval [0,1] 


(X n x' m - X m X' n y = (A m - A n )A m A, 


and using the boundary conditions. The appearance of the scalar product (III. 151) suggests 
to look for Hilbert spaces adapted to it. We consider this issue in the next section. 


We want to point out here that, at this stage, the position of the point particles are the 
limits X(j) = ]irn x ^j A (x). However, as will be shown in the following, there are technical 
reasons to introduce function spaces where generically this equality does not hold (instead, 
boundary conditions of the type (IA.8|) are satisfied). 



























II.3. Alternative Lagrangian for the system 


In order to get a precise Hamiltonian formulation for the dynamics of the system at hand, 
properly identify and characterize its degrees of freedom and deal with the delicate analytic 
and geometric issues posed by the presence of boundaries, it is most appropriate to use the 
GNH geometric algorithm developed in [l3-16|. A convenient starting point is to introduce 
a new Lagrangian defined on a manifold domain 24] of the tangent bundle of a configuration 
space, that we will take to be a real Hilbert space. In practice this means that we will have 
to extend the system somehow and also consider held configurations which are less smooth 
than the ones used in (III. 3D . We will require, nonetheless, that the solutions to (III. 1D - (III.6D 
are appropriately contained in those corresponding to the equations of motion of the new 
Lagrangian. 

In view of the results of previous section, instead of working with (III. 3D . it is natural 
to look for a generalized Lagrangian written in terms of the scalar product defined by a 
certain measure // and the associated RN derivatives. The hope -that will be realized- is 
that the equations of motion give rise to a standard Sturm-Liouville problem in terms of this 
derivative and also that the boundary conditions defining the elliptic operator that will play 
a central role in its solution are such that its self-adjointness (and other related properties 
such as the completeness of the set of eigenfunctions) can be readily asserted and proved. 

Let us consider then 


^ i r9r „ ,, _.du, . . , . _, du . d 2 u 

V lu€ L-\ 0,1] : 3— (x) Vx € [0,1], — is^-a.c., 


e Ll°. 1 ] 


and the Lagrangian L : T> X L^[0,1] —> M of the following form 


L(Q,v) = -<y,v) lt 


1/dQ d Q 

2 \ dp ’ dp 


-uj 2 (Q,Q)l 


(11.16) 


- £ <4 ( ( - 1)j |: W) - A v>w>) - \ £ 


16 ( 0 , 1 } 


16 ( 0 , 1 } 


where (■,•)/* is the scalar product with respect to the measure /i = ao^o + hi|(oi) + aiAi, 
denotes the associated RN derivative and /i-a.c. stands for “p-absolutely continuous”. 
In view of the scalar product (III.15D . one might naively expect that a.j = Hj, however oq, 
as well as A(i), have to be taken as non-trivial functions of the physical parameters of the 
model (see 21] and appendix [A] for more details, in particular equations (1A.8D and (1A.11D ). 

The equations of motion can be obtained by computing the first variation of the action 
derived from the Lagrangian. A straightforward computation gives 


SS(u) 



(u, du)^ + (A ^u, du - u; 2 (u, du) fl 

~ (( _1 ) Ja l7l(^ / ) - 7iOh>) ((-l) J ^(i) - Mj)u(jf) 

16(0,1} c 


where du G V so that the traces 7 j(du') and 7 j(du) are well defined (see Appendix IA.IID . 



























From this last expression we get the equations of motion 
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u — (A m — u> 2 )u = 0 x G [0,1], /i-a.e. (11.17) 

- AU)n(j) = 0 j € {0,1} . (11.18) 

ci/i 

They have the form of the 1+1 dimensional Klein-Gordon equation on the interval [0,1] 
subject to Robin boundary conditions written in terms of the RN derivative. It is important 
to remember that A M is not self-adjoint in D. However, the solutions to (III. 181) belong to 


V := < u G V : 


and in this domain A^ is indeed self adjoint (see section lA.IIID . 

In order to see that these equations describe the same dynamics as (III. dll - (III. 61) we first 
notice that in the open interval (0,1) equation (III. 171) is simply the Klein-Gordon equation 


u 


u" + u 2 u = 0. If we write now the boundary conditions (III. 181) in the form given by 


(1A.8I) , plug the resulting expression into u(j) — A +u(j) + u 2 u(j) = 0, and use the relations 
that fix aij and A{j) in terms of the physical parameters for the problem we immediately 
obtain (]II.5I) - (III.6I) . This completes the proof of Proposition Q] 


II.4. Hamiltonian formulation 

The equations (III. 171) and (III. 18[) derived from the Lagrangian (III.16D _can be understood 
as a particular case of the abstract wave equation (see discussion in [l3j]). By using the 
results described there, it is possible to directly get both the expression of the Hamiltonian 
vector field and the manifold domains where the dynamics takes place. 

Proposition 2. The Hamiltonian dynamics takes place in the second c/as.*!! ( generalized ) 
submanifold A4 2 = V x V of the weakly symplectic manifold (X4i = V x L^lo), where lu is 
the pullback to M. v of the strong, canonical, symplectic form on L ^ x L 2 r The dynamics is 
governed by the Hamiltonian vector field X = ( Xq,Xp ) : XI 2 —> Xii whose components are 

Xq(Q, P) = P (11.19) 

X P (Q,P) = -(u 2 -A,)Q, (11.20) 

and the Hamiltonian is given by 

H(Q, P) = 1<F, P)„ + 5 (f. f 1 + ^(Q, Q)» (11.21) 

+ Y. A(j)Q\j). 

16 ( 0 , 1 } ^ ^ 16 ( 0 , 1 } 


2 A submanifold IV A M of a presymplectic manifold is said to be second class if TA+ flTW = {0} 

where TN ± := {Z e TM\ N : w\ N (Z,X) = 0VX £ TN}, and TN_ := j*TN. 
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As we have mentioned before, this result follows directly from the abstract wave equation 


13], however, there are interesting technical details in the derivation of the Hamiltonian 


formulation by using the GNH algorithm. For completeness, we sketch here how this works 
for our system. 

We start by computing the fiber derivative defined by the Lagrangian (III. 161) 


FL-.VxLl^Ll* L 2 ; : (Q, V) (Q, (V, ■}„). (11.22) 

In order to conform with the standard notation we will write it as (Q , (P, ■) (U ) with P := V. 
By using now the Riesz representation theorem we can simply consider, as in the standard 
case of the scalar held with Dirichlet or Robin boundary conditions [l3j, that the fiber 
derivative is P x L 2 l —> L 2 x L 2 with FL(Q , V) = (Q, V ) and the primary constraint 
manifold is Adi := T> x L 2 taken as a generalized submanifold of L 2 x L 2 . 

The space L 2 x L 2 carries a canonical strongly nondegenerate symplectic form (inherited 
from the cotangent bundle L 2 x L 2 *) given by 

fi (Q,P)((9i>Pi)> (Q2,P2)) = (9i,P2) m - {q2,Pi)n (H-23) 


where we have Q,P,qi,pi G L 2 . The pull-back of fl to Adi is the weakly symplectic form 
cu := FL*Vl given by 


U(Q,P){(quPi), (q2,P2)) = (gi,P2)» - (Q2 ,Pi)v , (11.24) 

where we have now Q, qi e T> and P, pi e L 2 . It is interesting to mention that the “boundary 
terms” of the scalar product give rise here to “boundary terms” in the symplectic form. 
However we will see that this does not imply the existence of boundary degrees of freedom. 
In order to obtain the Hamiltonian on Adi we compute the energy 


HoFL{Q,V) = {V,V) li -L{Q,V). 


(11.25) 


This expression fixes the values of the Hamiltonian F[ only on the primary constraint sub¬ 
manifold Adi, however, this is the only information that we need to proceed with the GNH 
algorithm. From (III.251) we find that H : Adi —> M is given by Equation III. 211 

On the primary constraint submanifold Adi, vector fields are maps X : Adi —> Adi x Adi : 
(Q, P ) ^ ((Q, P), (X q {Q, P),X p (Q, P)), such that X q {Q, P) e V and X P {Q, P ) G L 2 . We 
have then 

(ixu)(Q,P)(q,p) = (XqiP)h ~ (q,X P )n. 


We must hnd now a submanifold Ad 2 and an injective immersion Ad 2 A Adi that allows 
us to solve the equation 

(hw - d H) | J2 (m 2 ) = 0 . (11.26) 

Notice that this will require us to identify Ad 2 as a subspace of Adi and also to specify its 
topology (by giving, for instance, a scalar product on it). 

The resolution of Equation (III.261) is relatively long but straightforward (see [l3[ for 
several similar computations). The final result is that Ad 2 = T> xP and the Hamiltonian 
vector held is given by 


Xq(Q,P) = p 
Xp(Q,P) = -(£ 2 -A M )Q. 


(11.27) 

(11.28) 
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The injection j 2 : M 2 — > M\ is just given by the inclusion map and is continuous if we 
use the natural topologies defined by the scalar products in M\ and M. 2 . The Hamiltonian 
vector field ( Xq,X p ) is also continuous in the same topology. 

By using the same kind of argument employed in [l3j, it is straightforward to show that 
the closure of T> x T> in T> x L 2 satisfies clp xi 2 {T> xD) = X? x L 2 . We then conclude that the 
Hamiltonian vector field X is tangent to Xi 2 so the GNH algorithm stops. We finish this 
section with a comment. Although it is not obvious at first sight, the Hamiltonian (111.211) 
is positive definite when restricted to M 2 . 


III. FOCK QUANTIZATION 

The fact that the positions of the point particles attached at the ends of the string are 
not independent physical degrees of freedom actually suggests that the Fock space for this 
system will not have the form of a tensor product of different Hilbert spaces associated with 
the masses and the string. Although this is quite natural from this perspective there are, 
however, physical questions that come to mind. For instance, if we think about this model 
as two masses connected by a physical spring (with “internal degrees of freedom”), how does 
one recover the situation where the string just models an ideal spring? What is the origin 
of the L 2 (M) ® L 2 (M) Hilbert space that one would use to describe this system? As a first 
step towards addressing these questions it is important to understand in detail why the Fock 
space does not factorize. We do this in the following. 


III.l. Construction of the Fock space 

An accepted way to quantize linear systems relics on the construction of a Fock space. 
In the specific case of quantum field theory in curved spacetimes the relevant details can be 
found in [ 13 ]. The starting point is the real vector space S of the classical (smooth) solutions 
to the field equations. By introducing a complex structure (among the many available) a 
complexified version of S is defined and endowed with a sesquilincar form obtained with the 
help of the symplectic structure. Taking a subspace of “positive frequency” solutions, the 
sesquilinear form defines a proper scalar product (•,•)+• By completing the complex vector 
space of positive frequency solutions w.r.t. this scalar product we obtain the one particle 
Hilbert space f). Finally, the Hilbert space of states of the quantum field theory is the Fock 
space T-L := .F(f)). 

In our case, the Hamiltonian description of the preceding section has produced a linear 
manifold domain of L 2 x L 2 given by Xi 2 = T>x P, where the classical Hamiltonian dynamics 
takes place, and a Hamiltonian vector field X tangent to the closure of M 2 . Notice that we 
can pullback the canonical symplectic form in phase space 

ty<?,P)((9i>Pi)> (Q2,P2)) = (qi,P2)n ~ (q2,Pi)» (HI.1) 

to M 2 (where (Q, P), (q^Pi) G Xi 2 ). We introduce now a complexification of this vector 
space M 2 and use the symplectic form to define a scalar product. Vector addition is defined 
componentwise as the standard sum of real functions and multiplication by complex scalars 
is defined by introducing the complex structure: 

J : M .2 x M 2 —X M 2 x M 2 '■ ( Q 2 ,P 2 )) H > ((—Q 2 , —P 2 ), (Qi,Pi)) 


(III.2) 
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and requiring that 


{a + bi) ■ V := (aI + bJ)(V) 


(HI 3) 


for a,b gR and V G XI 2 x M. 2 - In this case one can think of the elements in the complexified 
vector space as complex functions in XI 2 with the standard sum and multiplication by 
complex scalars. The complexified symplectic form is the straightforward extension by 
complex linearity of (IlIl.lD 


((<?!,Pi), (?2,P2)) := (qi,P2)p - (?2 ,Pi)m- 


The integral curves of the Hamiltonian vector held X can be identified with the space of 
solutions to the equations of motion. As X is defined in terms of the elliptic operator A M , 
its eigenvalues X n and normalized eigenfunctions Y n satisfying A ^Y n = X n Y n (see appendix 
IA.IVI) . play a relevant role in the following. It is important to point out here that these are 
not necessarily classical, regular solutions (for instance, they are not smooth with compact 
support). Furthermore, as the evolution is a symplectic transformation, the symplectic form 
can be pulled back to this space in the obvious way. In particular, the complexified solution 
defined by the (complex) Cauchy data (Q,P) (u(-, 0), u(-, 0)) G Adf at t = 0 can be 

written in the form 


u(x, t ) = 


\ 'y ^ f e ity/w 2 -x n 

n ' 




e ~ity/uj‘ 2 -\ n 



with Q n = (YnjQ)^ G C and P n = (Y^P)^ G C. Notice that, as mentioned in section UL2l 
our assumption that u > 0 implies that X n < uj 2 . In order to select a subspace of positive 
frequency solutions we require that 


iP + yj 2 - <3 = 0 . 


(111*4) 


This condition is equivalent to 


Qr 


Pr, 


x/uP 


Xr 


= 0, Vn. 


It is now straightforward to see that, when equation (1III.4D holds, we can indeed define a 
scalar product in this “positive frequency part” of the solution space as 


(Q (1) , Q (2) }+ := -AI C ((Q« p (1) ), (Q (2) , p (2) )) = 2(gw, ;q (2) ), . (III. 5) 

At this point, the only remaining step to finish the construction of 1) is to Cauchy complete 
in the norm defined by the scalar product. In terms of the Fourier coefficients of Q ^ and 
Q^ in the orthonormal basis that we are using the preceding scalar product becomes 

<Q (1) , Q&) + = 2 Vu 2 - A n (Y^QW)„(Y n , Q^),. 

n 

Notice that the L^-orthonormal basis {Y n } leads to an orthonormal basis of t) consisting of 


Zn ' (Y^YP 1 / 2 ^ V2(^~x n yp 


Y n = 


Y n 
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hence, the one particle Hilbert space can be identified with 


= 


i ^ ^ 

\ n 


Z n : ij) n e C, \^Pn | 2 < °o 


Finally the Hilbert space H of our quantum field theory is given by the symmetric Fock 
space over 1), "H = F( f)). 


The standard procedure described in [17|] can be used to introduce creation and annihila¬ 
tion operators, quantum fields and the quantum Hamiltonian H that generates the quantum 
dynamics of the system. In particular H is given by the lift to the Fock space of \Juj 2 — A 


A 4 ' 


Notice that the unitarity of the quantum evolution is guaranteed by the self-adjointness of 
the operator A„ in T>. 


III.2. Factorization of the Fock space 


In this section we will discuss the impossibility of factorizing H = H, n ass w © String in a 
natural way, i.e. with “factors” associated, respectively, with the point masses and the string. 
For the sake of the argument let us see what would happen if the Hilbert space for our model 
were the Fock space over l?. As discussed in appendix IA.II the space L 2 is isomorphic to 
R © L 2 (0,1) © R where each R is associated with a point mass. Using the properties of the 
Fock construction it is straightforward to see that F(Lj t ) = J-'(R) © J r (L 2 (0,1))©.F(R) and, 
hence, if this were the case, we would have a Hilbert space corresponding to each of the 
point masses contributing a factor to the Fock space. 

However we have to consider 1), instead of L 2 , and we are prescribed to use the scalar 
product (■,■)_)_ involving the square root factor sju 2 — A^. This changes crucially the out¬ 
come and prevents us, in particular, to write b = C © String © C with the C subspaces 
associated with the boundaries as in the previous example. Indeed, if such a decomposition 
were available, the function 


F : [0,1] —>■ R : x (->• F(x) 


1 x = 0 

o x e (o, l] 


(III.6) 


considered in the appendix (see equation (IA. 12)1 and the discussion after it) would have to 
be in 1). However, it is straightforward to see that, in the limit n —> oo, the coefficients 
(F, Z n ) + satisfy 

(F, Z„> + = ddg_ + o (n-^) , 




so they are not square summable and, hence, F ^ fy 

The impossibility of achieving the previous decomposition immediately implies the non- 
factorizability of the Fock space built from f). Of course, this does not exclude the possibility 
of finding such a factorization in other ways, for example for every decomposition of the type 
f) = ()i © f )2 we would have 


n = J-(fj) = F(bi © h) = -F(f)i) <8> F %), 

however, they are definitely not obvious from the present perspective and, in particular, 
there is no reason a priori to associate the factors to the point masses. 
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IV. CONCLUSIONS 


We have studied the Hamiltonian formulation and Fock quantization for a 1+1 dimen¬ 
sional model containing both fields and point masses. The combined description of different 
types of physical objects poses some interesting and non-trivial questions related to the 
proper characterization of the physical degrees of freedom both at the classical and quan¬ 
tum levels. Some natural questions for these systems are: Is it possible to talk about 
independent degrees of freedom associated with the masses and the fields? or, is it possible 
to split the quantum Hilbert space of the system in the form % = ® 'Hfieid? One 

of the main results of the paper is the proof that the Fock quantization of such compound 
model leads to a Fock space that cannot be written in a natural way as the tensor product 
of Hilbert spaces associated with the masses and the field respectively. We want to em¬ 
phasize that this result comes about after a careful discussion of the spaces of solutions to 
the equations of motion and the details of the construction of the Fock space, in particular, 
the scalar product appearing in the complexified solution space. This lack of factorization 
must be understood because one of the assumptions used in the description of compound 
quantum systems is that their state spaces are tensor products of Hilbert spaces associated 
with each subsystem. This statement is singled out by some authors as the zeroth postulate 


of Quantum Mechanics 25), [26|. Of course, in the present case our construction does not 


violate this postulate but, rather, shows that identifying subsystems in a proper way is not 
straightforward, even at the classical level. 


From a technical point of view the method that we have employed relies on the intro¬ 
duction of Hilbert spaces endowed with scalar products defined with the help of modified 
measures describing the coupling of the fields with lower dimensional objects (point masses 
in the present case). This is useful both at the classical and quantum levels because the 
identification of the relevant functional spaces is simplified in a significant way. One of 
the key elements of our approach relies on the ideas developed by Evans [2lj] to deal with 
modified Sturm-Liouville problems of the type considered here. It is especially important 
to work with elliptic, self-adjoint operators (generalized Laplacians) with spectrum and as¬ 
sociated eigenfunctions that can be mapped to the ones appearing in the resolution of the 
non-standard eigenvalue problem that determines the normal modes of the system. These 
generalized Laplacians play a fundamental role in the Hamiltonian formulation of the model 
and its subsequent Fock quantization. 


The Hamiltonian formulation that we have developed starts from a Lagrangian written 
in terms of the natural objects: measures, the scalar product defined with their help and the 
associated Radon-Nikodym derivatives. A striking feature of the approach that we follow 
is that although the constants appearing in the scalar product are functionally dependent 
on the physical parameters they are, generically, non-trivial functions of them. It is also 
important to emphasize the fact that the appropriate Lagrangians are non-trivial when 
written in terms of these objects because of the necessity to have appropriate self-adjoint 
operators (in particular it is very important that the normal modes of the system can be 
interpreted as eigenfunctions of the self-adjoint operator A M with eigenvalues given by the 
normal frequencies). 

We have described the construction of the Fock space by relying of the methods custom¬ 
arily used in the context of quantum field theory in curved spacetimes (in particular those 
described by Wald in (T3] ) - An advantage of combining those methods with the Hamiltonian 
formulation that we are using here is the possibility of having a precise and explicit char- 
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acterization of the solution space for the equations describing the dynamics. Basically the 
only difference with the present case is due to the fact that the relevant elliptic self-adjoint 
operator is not a Laplace-Beltrami operator but one defined with the help of a measure with 
singular contributions at the boundaries of the region where the fields are defined. 

We have illustrated the ideas developed in the paper with a simple 1+1 dimensional 
model but they can be extended without difficulty to other more complicated systems, in 
particular, our approach can be used to deal with higher-dimensional models and is flexible 
enough to allow the coupling of different types of low-dimensional objects not restricted to 
point masses. The model that we have studied provides an interesting way to interpolate 
between different types of boundary conditions (by, for instance, considering the limits of 
small or large masses at the boundary). It is important to highlight here that, the present 
methods can be of use not only for linear systems but also for other non-linear ones that 
are obtained by consistently adding interactions to them (for instance, gauge theories). By 
working with fields defined in unbounded space-time regions these techniques can also be 
adapted to the study of quantum dissipative systems (in the spirit of the Caldeira-Legget 
or Rubin models) and decoherence. 

A final comment concerns the interpretation of the current system in terms of so-called 
particle detectors, which provide a clear-cut definition of a particle in general situations, 
such as curved spacetimes and non-inertial frames. In this approach, the particle detector 
is a system with additional degrees of freedom that interacts with the field according to 
some specific interaction Hamiltonian, as in e.g. The effect of the interaction is read 

in the evolution of the states of the Hilbert space, which is a tensor product, and in which 
excitations and de-excitations in the sector of the detector according to its interaction with 
the field are interpreted as particles. For an in-depth review of the state of the art see 27 


As we have shown this type of factorization does not take place for such natural models as 
the one considered here so it is very important to understand and characterize the physical 
systems for which this is possible and the ensuing implications for the measurement problem 
in quantum mechanics. 


Appendix A: Hilbert spaces, modified measures and Radon-Nikodym derivatives 

A.I. The Hilbert space L 2 


The implementation of the GNH algorithm requires the introduction of appropriate 
Hilbert spaces. After realizing that a scalar product of the type (III.1511 plays a natural 
role in the eigenvalue problems that crop up in the study of the system that we are consider¬ 
ing here, it is natural to construct Hilbert spaces endowed with it. To see that this is not a 
completely trivial task, it suffices to realize that (11 1 .1 51) does not make sense for elements of 
L 2 ( 0 , 1 ) as they do not have well defined boundary values at x — 0,1 (because points have 
zero measure with respect to the standard Lebesgue measure). A possible way out could be 
to consider other types of Hilbert spaces, for instance the Sobolev space H l { 0,1). The extra 
regularity of the elements in Sobolev spaces makes it possible to define their boundary values 
by using the so called trace operator 7 . The elements of the H l (I ) Sobolev space defined on 


an interval / Ci have representatives given by continuous functions on its closure / 28 


Their values at the boundary (or their limits for that matter) are the traces that we denote 
as 70 , 71 because the boundary of an interval is disconnected. 

In the present context it is better to avoid demanding so much regularity so we consider, 
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to begin with, spaces with elements satisfying minimal regularity requirements. The simplest 
example of such a space would be L 2 [ 0 , 1 ] := 1 x L 2 ( 0 , 1 ) xK consisting of elements that 
we will denote as v := (v 0 , u, tq). This is a Hilbert space with the scalar product 

( v , w)l := a Q v 0 Wo + otiViWi + / v ■ w c\/i L (A.l) 

J[ o.i] 

provided that ay > 0. Notice that the completeness of L 2 [0,1] is a direct consequence of 
that of L 2 ( 0 , 1 ) and M. 

A useful alternative way to describe this space makes use of a suitable measure p. Let 
us consider the Borel u-algebra 03([0,1]) and the measure p = ao^o + /Gl|(o,i) + cqhi, where 
the measure Sj is defined for any A e 03([0,1]) by 

6 ^ = {lV 3 l A A < A - 2 > 

Let us consider now the real Hilbert space of square integrablc functions 

L 2 [ 0 , 1 ] := i / : [ 0 , 1 ] —> M : / is p-measurable , / / 2 dp < +oo 

l J[o,i] 

(as usual, functions equal p-a.e. are identihed) endowed with the scalar product 

(u, v) M = / u ■ v d/i . (A.3) 

J[o,i] 

It is straightforward to see that the map 

4 : LJ[0,1] —> L 2 [0,1] 

/ (/(0),/!(„,!),/(I)) 

is a Hilbert space isomorphism, hence, in the following we will view elements of these Hilbert 
spaces as square integrable functions (w.r.t. the measure /i) defined on the closed interval 
[ 0 , 1 ] or as elements of Ixi 2 (0, l)xl with the scalar products defined above. Whenever 
no confusion may arise we will use the shorthand L 2 to denote L 2 [0,1]. 



A.II. Absolutely continuous functions and Radon-Nikodym derivatives 

A useful way to write the Lagrangians that we use in the paper makes use of the Radon- 
Nikodym (RN) derivatives of appropriately defined functions. They also play a central role in 
the description of the relevant elliptic operators that appear in the Hamiltonian formulation 
and the construction of the Fock space. Here we briefly review the central concepts and give 
a list of properties that are used throughout the paper. 

Given two measures p and v over f8([0,1]), we say that v is p-absolutely continuous (/i- 
a.c. usually denoted as v <C p) if whenever p(A) = 0 for some A e 93([0,1]), then v(A) = 0. 
The RN theorem states that over finite measure spaces, such definition is equivalent to “weak 
p-differenciability”, which means that there exists a /j-measurable function / e L 1 (/i) such 
that 

”f( A ) = f / d F- 

J A 
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The function / (usually denoted by f — (A) is known as the RN derivative and it is unique 
in the sense that any other function that satisfies the preceding properties is equal to / 

/i-cl. 6. 

If we want to /^-differentiate a function F : [0,1] — > R w.r.t. p = a 0 ^o + / 7 l |( 0 i) +«i^i we 
have to associate to F a /i-a.c. Lebesgue-Stieltjes measure //p. To this end we define 

i'f = (f( 0+) - F(0)),So + F'iil | (01) + (f( 1) - F(l-))ii 

where F^^ is a - c - i n the usual calculus sense (so that F'/x^|( 01 ) is pl- a.c.). Under these 
conditions F is differentiable /x^-a.e. and has well defined limits F(0+) and F(l —), not 
necessarily equal to its values at the boundary F(0) and F(l). In fact, considering that the 
a.c. functions over (0,1) can be seen as elements of the Sobolev space FT 1 (0,1) (functions 
of L 2 ( 0 , 1 ) with distributional derivative in L 2 ( 0 , 1 )), and also that no restriction over the 
boundary { 0 , 1 } is placed by the condition of being /i-a.c. continuous, we conclude that the 
vector space of p- a.c. functions H^[ 0,1] is isomorphic to R x H 1 { 0, l)xR (compare with 
L 2 [0, l]=lx L 2 ( 0,1) x R). An element F G H^[ 0,1] has boundary values F(0), F( 1) and 
traces 70 (F) = F(0+), 71 (F) = F(l—) which are generically different from F(0) and F(l). 
This is related to the fact that, in general, aq ^ p 3 . In our model the boundary values 
of functions (such as F(0), F(l)) play a secondary role from the physical point of view as 
the positions of the point masses will be defined by the traces, but they play a central role 
from the mathematical point of view because they are crucial to prove that the fundamental 
Laplace operator of the model is self-adjoint. 

Finally we define the RN derivative of F G H^[ 0,1] as := Notice that if hh is in 
Hjf 0,1], we can apply again the RN theorem to obtain another derivative that we denote 

777 := TI (ti)' The ex Phcit f° rm °f the RN derivative of F G 0,1] can be immediately 
obtained from the previous result by comparing the expressions of up and p 






(x) = F'(x ), pl- a-e. in (0,1) which is just the RN derivative of F w.r.t. pp- 
dp 


It is important to notice that the standard Leibniz rule does not hold, in general, for RN 
derivatives though a modified rule exists. Indeed, by defining K : [0,1] —> R with K | ^ ^ = 0 
and K(j ) = (—l) J 'aq for j G {0,1} we have: 


d(FG) 

dp 


(x) 


dF 

dp 


dC 1 H W dC 1 

(x)G(x) + F(x)^-(x) + K(x) — (x)-^(x ), p- a.e. in [0,1]. (A.4) 

dp dp dp 


This expression will be very useful whenever we have to perform integrations by parts of 
RN derivatives. 
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A.III. The fundamental Laplace operator 


Next we will introduce a generalized Laplace operator, defined in terms of the RN deriva¬ 
tive, that plays a central role in the paper. Let us write 


V : = { u e L 2 [0,1] : 3^(z) Mx G [0,1] , ^ isp-a.c., ^ G L 2 [ 0,1] 


dp 


dp 


dp 2 


; d U 


V ■= jw G V : = 0 

A M : C I^[0,1] -A L^[0,1] : « >->• (1 + C) 


d 2 w 

dp 2 


(A.5) 

(A. 6 ) 

(A.7) 


where C : [0,1] —* M and A(j) G M for j G {0,1} are to be determined together with aq. 
One condition that we have to impose is that we should recover the equations of motion, 
in particular for the lateral limits at the boundaries. Fulfilling this requirement does not 
completely fix these constants, but puts the eigenfunctions of A M in correspondence with the 
ones of the original problem. Hence, if we make A M self-adjoint, then the set of eigenfunctions 
will form a complete set of orthogonal functions w.r.t. the scalar product (v) M . 

So let us compare equations (III. lOP - OII.12D with the eigenvalue problem A = A Y 


(defined in a subspace of L 2 ). First, (III. 101) requires that Cj( 01 ) = 0 because then A a 


is the usual Laplacian (second derivative) when restricted to the interval (0,1). Taking 
C(j) = A(j)otj makes A M symmetric (i.e. (A — (u, A^w)^ = 0 ) provided we restrict 
its domain to T>. In fact using (1A.4(1 we have: 




E 


d 2 u (, , w d w ... 

I (— lyoij— -G(j)ra 


dp 


= - /- -\ 4 . y, mu(i)wu) 


\ dp ’ dp 


P J 6 { 0 , 1 } 


where we have used that the boundary conditions (—l) J ^)(j) — A(j)u(j) = 0 and our 
particular choice of C(j), leading to a symmetric expression in u and w. 

It is important to notice that ussing the expression given for the RN derivative at the 
boundary, we can show that the Robin-like boundary conditions (—l) 2 A)(j) — A(j)u(j) = 0 
are equivalent to 

7 j(u) = (1 + ajA{j))u{j) . (A. 8 ) 


fij Ao;? 


Now taking A (j ) = —AA and using the previous equation for dL and for 7 j(F), we 




d [a 


obtain that A = XY on ( 0 , 1 ) is equivalent to: 
Y" = A Y , 


lA Y> ) = ( _1 ) J Vji 1 ~ ajUjAujj) X + nj Aw. W j(Y) ■ 


(A.9) 

(A. 10 ) 


If we finally require aq(l — aqpjAcu 2 ) 2 = pj, equations (1A.10I) become equations (III.llfl - 
(111.12(1 . recovering thus the original problem. Clearly from this condition and the fact that 
fjbj > 0 , we see that all the possible solutions for oq are positive, so they can be used to 
define scalar products of the form flA.ll) . (IA.3l) . A very simple case corresponds to Aujj = 0 
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(some sort of resonance that cancels out the effect of the springs and the string) for which 
a.j = /ij and A(j) — 0. It is important to realize that, as we have already mentioned before, 
the positions of the point masses are given by the traces 7 'j(Y) (which can be identified with 
the X(j) of section Hi.21) . however, the values Y(j) that appear in the Robin-like boundary 
conditions (1A.8h do not have a direct physical interpretation as far as we know. 

At this point, we have only managed to obtain a symmetric operator A^ whose eigen¬ 
functions are in correspondence with the ones of the original problem, but in fact, it can be 
shown that the operator A M is self-adjoint in T> by using the method described in 21 ]. As 


mentioned before, this implies that the set of its eigenfunctions can be turned into a Hilbert 
basis of the Hilbert space T> showing in particular, that it is complete. A remark is in order 
here, if u,v G V we have, using equation (1A.8I) . that: 


(u,v)n= Y a Md) v (j) + («d') = Y /HTj'MTjO) + (u,v) = i(u,v )). (A.11) 

16 ( 0 , 1 } 16 ( 0 , 1 } 

This establishes a connection between the results in the present section with the argument 
leading to (III. 151) (notice that the boundary values of the eigenfunctions appearing there can 
-and should- be interpreted as right or left limits or, in our notation as 7 ^(w), 7 ^( 1 ))). 


A.IV. Normal modes and Fourier coefficients 


An orthonormal (Hilbert) basis of L 2 [ 0,1] can be constructed by using the eigenfunctions 


of the self-adjoint operator A^ defined in (1A.7|) in the domain T> given in (IA. 6 |) . These 
eigenfunctions, in turn, can be computed from the solutions to the eigenvalue problem 
considered in (III. 1QD - (1II. 121) by defining their values at x G {0,1} in such a way that the 
Robin-like boundary conditions, or equivalently equation (lA.Sp . are satisfied. The elements 
of the Hilbert basis will, hence, be of the form {Y n } (in one to one correspondence with 
the eigenfunctions {X n } introduced on section Hh2l) with Y n h 01 ) = X n /g n (in particular 
lj(Y n ) = X n (j)/g n ) and Y n (j) = (1 - ajPjAu>f)X n (j)/g n . The factor g n guarantees that Y n 
is normalized w.r.t. the scalar product (v) M . It can be explicitly computed to be 


2 ^ / a ~2 I /I I \ 2 1 2/2 a ~2\2 1 2/ 2 1 a ~2\ U n A ho( W n 

9 n ~ 9 ( h 0 AWq + (1 + h 0 )A + ho \ U n ~ ^o) + hlhAi + ^ 1 ) 2 i— a~ 2\9 

2 V Wn + hRWn-AwJ- 


The asymptotic behavior of 1 /g n when n —> cxd can be obtained by using equation (III. 141) : 


-i n /2 1 i n t -4} 

9n =—^TA + 0 \ n ) 

Pq ir z n z v ' 


An element F G L 2 [ 0,1] can be expanded as usual as F — ^2(Y n , F)JY n . It is interesting 
to notice at this point that functions / : [ 0 , 1 ] —» R supported at the boundary { 0 , 1 } are 
non trivial elements of L^[ 0 , 1 ] and hence can be expanded in the basis introduced above. 
In the particular case of the function F G L 2 [ 0,1] given by 


F : [0,1] — y M : x hg F(x) = 


x = 0 
x G (0,1] 


(A. 12 ) 
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we have (Y n ,F) M = a 0 Y n ( 0) = a 0 (l - a O R O Auj%)X n (0)/g n = ^n 0 a 0 X n (0)/g n . In the dis¬ 
cussion of the Fock quantization of the model we need the asymptotic behavior of these 
coefficients for large values of n. This can be easily obtained by using the fact that in this 
limit only the eigenvectors associated with negative eigenvalues matter. The asymptotic 
behavior of g~ l and the one of u> n given in equation (III. 141) lead to 


(Y n ,F), 


V To a o 


* n (0) 

9n 


\/ To a o 


^n 

9n 


a/2/A 


0«0 


2- + O (n~ 2 ) . 
nn 
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